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We study an electrically-controlled quantum spin Hall antidot embedded in a two-dimensional topo- 
logical insulating bar. Helical edge states around the antidot and along the edges of the bar are 
tunnel-coupled. The close connection between spin and chirality, typical of helical systems, allows 
to generate a spin-polarized current flowing across the bar. This current is studied as a function of 
the external voltages, by varying the asymmetry between the barriers. For asymmetric setups, a 
switching behavior of the spin current is observed as the bias is increased, both in the absence and in 
the presence of electron interactions. This device allows to generate and control the spin-polarized 
current by simple electrical means. 

PACS numbers: 71.10.Pm, 72.25 -b, 73.23.-b 



I. INTRODUCTION 

Two dimensional topological insulators showing 
quantum spin Hall (QSH) effect were recently sub- 
jected to intensive theoretical and experimental research. 
The theoretical prediction by Bernevig et a/P, about the 
possibility to realize the QSH effect in Hg-Te quantum 
wells, was experimentally proved by Konig et al.^ ^ who 
showed that these systems are characterized by gapless 
edge states at zero magnetic field. The absence of time 
reversal-breaking mechanism ensures the robustness of 
these edge states from backscattering.^ ^ The main prop- 
erty of QSH effect is that the edge states are helical,^ with 
locking between spin and momentum degrees of freedom. 
Several theoretical proposal were made in order to con- 
firm the connection between spin and chirality. In partic- 
ular, quantum point contacts (QPCs),^ extended con- 
tacts, inter ferometricsetups,IISHiHl ^s well as scanning- 
tunneling microscopy^^^^" turned out to be appropriated 
for this task. Recently, quantum dot realized in Hg-Te 
quantum wells have been theoretically studied. ^^^^^^ The 
fast development in the experimental techniques will soon 
allow to verify these theoretical proposals. 
The important role played by the spin degrees of freedom 
makes these system fascinating for future applications. 
The ability to generate and control spin currents would 
open the possibility to realize promising devices for the 
development of spintronics.^^ Recent proposals of devices 
for spin manipulaton, based on two-dimensional topo- 
logical insulators, rel y on A haronov-Bohm and Fabry 
Perot interferometers ,1^1^'^^'^ the application of a mag- 
netic field at a pn junction^^ and tuna ble wavefunction 
overlapping in nanoconstrictions. ^^'^^ ' However, simple 
two-terminal configurations are often not able to induce a 
global spin-polarized current flowing through the sy stem, 
unless magnetic fields or ferromagnets are present .1^^^ 



The aim of this paper is to analyze the possibility of gen- 
erating and controlling spin-polarized current in a two- 
terminal QSH antidot system. This geometry was al- 
ready proposed for detecting peculiar properties of frac- 
tional charges in fraction quantum Hall bar.^^ The in- 
terest related to this device is that the spin-polarized 
current can be controlled only by electrical means, thus 
overcoming the need for magnetic materials. 
The antidot is embedded in a QSH bar and coupled, via 
tunnel barriers, to the edges of the bar. We will show 
that this setup, by virtue of the close link between spin 
and chirality, allows not only to generate, but also to con- 
trol the spin-polarized current along the edges of the bar 
in a very simple and efficient way, both in the absence 
and in the presence of electron interactions. We will show 
that the asymmetry between the tunnel barriers is a key 
ingredient in order to produce spin-polarized current. 
The paper is organized as follows. In Sec. |TI|we introduce 
the theoretical model for the device. In Sec. Hill we de- 
scribe the master equation approach adopted in order to 
calculate the transport properties and we define the rela- 
tions between the current through the antidot and along 
the edges of the QSH bar. Section [IV| contains the main 
results of our research. The spin-polarized current is an- 
alyzed as a function of the external voltages, by varying 
the asymmetry between the barriers and the strength of 
the electron interactions. Section |Vl is devoted to the 
conclusions. 



II. MODEL 

The system is schematically shown in Fig. [l] An anti- 
dot is created inside a QSH bar in a two terminal setup. 
Gapless helical edge states appear along the edges of the 
bar and around the antidot. We model the coupling be- 
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Figure 1. (Color online) Schematic representation of an anti- 
dot embedded in a QSH bar. The edges of the bar are biased 
in a two-terminal configuration. Helical edge states appear at 
the boundaries between the QSH phase and the vacuum. Full 
red (dashed blue) lines represent spin up (down) electrons. 



tween the antidot {ad) and the upper (A) and lower (B) 
edges by tunneling barriers. The system is driven out of 
equilibrium by applying a bias potential V. The chemi- 
cal potential in the antidot region can be tuned through 
a gate voltage Vg. The full Hamiltonian of the system 



H = Ha + Hb + Had + 



(A) 



H 



Hv 



(1) 



consists of the sum of the leads {Ha + Hb) and the an- 
tidot (Had) uncoupled Hamiltonians, the tunneling bar- 
riers {h[^^ + ) and the coupling of the system with 
the external circuit {Hy). 

The helical edge states around the antidot are de- 
scribed in terms of helical Luttinger liquid (HLL)^ {h = 

1) 



Ha 



ad = -iVF I dx (t/'L^^^xV'ad.t " V-L.i^x^ad.i) 



(2) 



where vp is the Fermi velocity, L is the circumference 
of the antidot and g is proportional to the screened 
Coulomb repulsion. The bosonization technique^^ allows 
to write the electron operator 



"^adA^) 



(3) 



in terms of the bosonic field (/>cr(^), where the sign + 
(— ) is for spin up a =t (spin down a =1) electrons. 
The chiral index is omitted, since it is directly linked 
to the spin. In Eq. (|3|, represents a Klein factor, 
is the excess number of electrons with respect to a 
ground state number no, a, « is a short length cut-off and 



kp = 7r(no,^ + no^i)/L is the Fermi momentum. The 
Hamiltonian ([2| can be diagonalized and reads 



u=l,2 q>0 



En 

2 



(n 



-,2 , Eg 2 



(4) 

Here, dq^jy are annihilation operators of the collective 
plasmon excitations, while n = -\- and s = — 
represent the total excess charge and spin (a part for di- 
mensional factors — e and h/2) with respect to the ground 
state value no and sq. The capacitive coupling between 
the antidot and the gate is taken into account through 
the term in Ug = VgCg/e.^^ Due to boundary condi- 
tions, the wave number is quantized, q = 27:my/L and 
^v{(l) = ^z/Co/i^, with eo = 2t:vf/L the level spacing 
and rriy positive integers. Three different energies char- 
acterize the plasmon, the charge and the spin excitations 
respectively 
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The coefficient K = {1 ^ g / {iivf))~'^ encodes the 
strength of the interactions. For the non interacting case 
(^ = 0,i^ = 1) one has e = 2En = 2Bs. Even in the 
absence of electron interactions, the contributions cor- 
responding to the charging energy {En) and to the spin 
addition energy {Eg) are finite, due to the discrete nature 
of the dot and the Pauli principle. In addition, despite 
our microscopic model provides quantitative estimates 
for these energies, several influences that occur in exper- 
imental setups are here neglected. Indeed, long range 
interaction effects and coupling with gates can lead to an 
effective charging energy greater than the one obtained 
in the microscopic model. More robust is the spin ad- 
dition energy, which is affected at most by the exchange 
part of the interactions. Therefore, we treat in the fol- 
lowing En as a free parameter with En ^ Eg. 
One of the peculiar features of QSH systems is the ab- 
sence of spin-charge separation even in presence of elec- 
tron interactions. Indeed, the plasmon excitations are 
characterized by a unique energy scale,'^ in contrast 
to t he us ual interacting spinful one dimensional sys- 

|-gj^g[33|35l 

Analogously, the Hamiltonian of the edges A = A, 5, 
supposed infinite, is 



H^ = 



J —oo 



where '0A,cr(^) is the electron field operator in the edge 
A. Despite the edges and the antidot can be subject to 
different screening mechanisms, in this work we assume 
a unique interaction parameter g for the electron inter- 
action. The application of a bias voltage V in the two 
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terminal setup is described by 
1" 



(SA - sb) , (7) 



where nA,cr is the number of extra charges on edge A and 
8\ = ^A,t~^A,i- The couphng between edges and antidot 
is described by pointhke tunnehng barriers 



h[^^ =txYl ^idJ^^)M^o) + h.c. (8) 



with xx^xo the points where tunnehng occurs and 
the tunnehng amphtude. The latter can be controlled 
through gate voltages placed above the barriers, thus al- 
lowing to increase or decrease their opacities. 



the thermal averages, denoted by the brackets (...), 
we assume that collective excitations are relaxed by 
extrinsic damping processes over a time-scale much 
shorter than the average electron dwell time through 
the system. As a consequence, the distribution {rrij^^i} 
is weighted with a Boltzmann thermal factor and a 
sum over all final collective plasmon states {mj^j} is 
then performed. A given process is completely specified 
by the initial quantum numbers of the antidot (n, s) 
and by their variation after the process (An,As) only. 

The tunneling rates can be recast in a more compact form 



where, provided that ksT <^ one haJ^ 



(13) 



III. TUNNELING DYNAMICS 

The transport properties are determined by the tun- 
neling processes between the edges of the bar and the 
antidot region^. The state of the antidot is specified 
by the excess charge n, excess spin s and distribution 
{rriiy} of the plasmons. The corresponding occupation 
probability is determined by the tunneling rate T^^^ 
from an initial state \i) = In^, s^, {m^^^i}) to a final one 
I/) = |n/, 5/, {mj,j}) across the barrier A. The sequen- 
tial nature of the tunneling processes restricts the pos- 
sible final states to n/ = + A^, Sf = 5^ + A^, with 
A^, As = ±1. At the lowest order in the tunneling ampli- 
tudes, the transition rates can be written in the compact 
form 



r 



(A) 



dte 



(A) 



f'Gad{t)gx{t), (9) 



where AU} 



(A) 



^^Aj,A.(^.^) = MA.,A.(n,5) T ^A,. (10) 



AU 



(A) 

A^,A, 



(n, s) with 



eV 



The sign ^ in Eq. (10) distinguishes tunneling across 
barrier A (— ) and B (+). The quantity /iA^,As (^^ ^5) is 
the chemical potential of the antidot without plasmon 
excitations 



1 A / 

- + An[n - ng 



MA^,A,(^,5) E, 

InEq. ^ we introduced 

Gjit) = (V'i.t(o,t)V'it(0'0)) = (V'i4(o,t)V'L(o,o) 



(11) 



(12) 

the electron Green's functions on the edges {j = A^B) 
and in the antidot {j = ad) computed for the system in 
thermal equilibrium at temperature T. In performing 



ap = e{p + 0+ f^'^^^\ peZ (14) 



e 2 



2tt 



r(C) \/3uj. 



C-i 



(15) 

with B the Euler Beta function, T the Euler Gamma 
function and 6 the Heaviside step function. In Eq. ([l3| 
we introduced the energy cut-off Uc = vp/iKa)^ the in- 
teraction dependent coefficients ^ = {K ^ 1/^) /2, the 
inverse temperature j3 = l/ksT and the characteristic 
rate 



.(A) 



I^aI^ 



1 



on the barrier A. We point to Ref. 
tailed derivation of Eq. (13). 



(16) 

I and [37| for a de- 



In the sequential regime, the transport through the an- 
tid ot can be approached by a master equation formal- 
^gj^| 36 | 38| r^Yie occupation probabilities of the states |n,5) 
of the antidot in the stationary regime satisfy 



E E 

.(A) 



r(A) 



(n,5)- 



Pn+A.,.+A.rLX,-A.(^ + An,S + A,)J =0 (17) 

together with the normalization condition ^ Pn^s = 1- 
Charge and spin tunneling currents through the antidot 
between upper and lower edges are 



jitun) ^ _g 



E 



A P 



,A^,As 



■ A^,As 



(n,5) 



■ A^,As 



jitun) ^ ^ 

n,s, An, Ag 



■ An, A 



Xn,s) 



■ An, A, 



(n,5) 
(18) 



(n,5) 
(19) 
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The peculiar helical nature of the system reflects in a 
close connection between these tunneling currents and 
the longitudinal currents flowing along the two edges of 
the bar. First of all, tunneling processes are responsible 
for a decrease in the total charge current flowing from 
left to right contact with respect to the uncoupled edge 
conflguration, i.e. / = — I^s- It has been shown"^^ 
that the backscattering current Ibs is directly related to 
the above spin tunneling current from edge to edge by 



/ (b) 



2e 



j{tun) 



(20) 



The spin-polarized current flowing along the edges of the 
QSH bar^^ Isp = | Scr=t,i i^^,^ ~ is zero in the 

uncoupled edge conflguration. However, when the anti- 
dot connects the two edges, a non-zero charge tunneling 
current is responsible for the generation of a longitudinal 
spin-polarized current 



T — — jitun) 



(21) 



Ordinary two-terminal tunneling geometries, realized by 
means of QPCs or extended contacts, are not able to gen- 
erate a spin-polarized current along the edges of the bar 
only by electrical means. ^ 

However, what we are going to show is that, under cer- 
tain conditions, the antidot acts as a spin filter on the 
longitudinal current, allowing to generate a spin polar- 
ization and to control it in a very simple and efficient 
way. 




Figure 2. (Color online) (a) Schematic representation of the 
antidot coupled to the edges of the bar through tunneling 
barriers, according to Fig. [l] The bias voltage V controls the 
chemical potential difference of the edges. The gate voltage 
Vg controls the chemical potential of the antidot. (b) Scheme 
of the states involved in the transport. Thick black lines de- 
note ground-to-ground state transitions, colored solid lines 
represent transitions involving zero-mode spin excited states, 
while dashed lines correspond to a transition involving a plas- 
mon excitations. The modulus of the initial (/) and final {F) 
spin states are denoted as |s/| \sf\ (spin states =bs are 
degenerate in our model). 



V ^ antidot excited states become also involved. 
In the symmetric case (77 = 1), from Eqs. ([lO|, ( pTj ) 

= Tf^_^ {n,-s). 
(n8|, pTl, implies 



and (13) follows that ^^^l /^^{'^•, s) 
This property, together with Eqs. 



that /, 



{tun) 



= = Isp at any V, so that longitudinal 



spin-polarized current along the edges cannot be gener- 
ated. 



IV. SPIN-POLARIZED CURRENT 

At low temperature {ksT <C co^eV)^ the transport 
occurs predominantly via spin up (down) electrons 
jumping into the antidot through barrier A (5), while 
spin down (up) electrons jump outside the antidot 
through barrier A {B). In Fig. [2|a) the two tunneling 
barriers are represented with different thicknesses. We 

(A) 

deffne here the unit of tunneling rate as Fq = Fq 
and introduce an asymmetry parameter r] so that 

Figure [2|b) depicts a scheme of different transport 
regions. Coulomb blockade regions appear at low bias 
where sequential tunneling is blocked due to energy 
conservation constraints. Tunneling becomes allowed at 
ffnite V when AU^^ ^ (n, s) < 0. Moving out of the 
Coulomb blockade region, the lowest-lying transitions 
are those involving the ground states of n and n + 1 
electrons (we assume here n even for deffniteness). 
The processes that load one extra electron into the 
antidot are |n,0) ^ |n + through barrier A and 

|n, 0) ^ |n + 1,-1) through barrier 5, while the ones 
unloading it are \n + 1,1) |^,0) through barrier B 
and |n + 1,-1) 1^,0) through barrier A. Those are 
represented as thick solid lines in Fig. |2|b). Increasing 



The situation with asymmetric tunnel barriers is 
however markedly different. In the rest of the paper we 
consider the case rj > 1^ the situation with r] < 1 being 
completely analogous swapping edge A with B and spin 
up with spin down. Figure [3] shows the spin-polarized 
current Isp for (a) non-interacting {K = 1) and (b) 
interacting {K = 0.57) electrons. Note that in the 
non-interacting case, collective excitations lines overlap 
with the ones involving zero-mode spin excited states. In 
the interacting case, transitions are smoothed, because 
the edges of the bar are no longer Fermi liquid-like. 
Furthermore, plasmon excitations are pushed to higher 
energies and, in general, do not overlap with the ones 
involving zero-mode spin excited states. Note that two 
transport regions can be identified, supporting either 
positive or negative spin-polarized current. 

Figure |4] shows Isp{V) for rig = 0.5, for (a) non- 
interacting {K = 1) and (b) interacting {K = 0.57) 
electrons. Differences in shape between the two pan- 
els are due to the effects of interactions discussed 
above. However, by varying bias, in both cases the 
spin-polarized current alternates between positive and 
negative values. The alternance in sign of Isp is accom- 
panied by a decrease in the average value of the spin of 
the antidot (s) = ^/'^'^n g sPn^s^ as shown in Fig. [4|c) 
and (d). As we are going to show, these two effects are 
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be captured by looking at the most relevant processes 
occurring in the system. Let us focus the attention on 
regions / and //, which are indicated in Fig. |2|a): in 
region / (//), the maximum spin state available is |s| = 1 
(l^l =2). Consider region / and the antidot in its ground 
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Figure 5. (Color online) Scheme of the most relevant trans- 
port processes occurring in (a) region / and (b) region // (for 
?7 > 1). See text for details. 



Figure 3. (Color online) Spin-polarized current Isp (units 
hVo 1 2) flowing along the edges of the bar as a function of V 
(units eo/e) and rig for K —\ and (b) K — 0.57. In both 
panels ry = 10, /csT/eo = 0.02eo, ^n/eo = 20, cJc/eo = 100 
and n = 0. 




Figure 4. (Color online) Plot of /sp (units of /iro/2, upper 
panels) and (s) (units of /i/2, lower panels) as a function of 

V (units of eo/e) for rig = 0.5, corresponding to the dashed 
yellow hne of Fig. js] with (a-c) i^^ = 1 and (b-d) K = 0.57. As 

V sweeps, /^p changes sign and |(s)| increases. In all panels 
T] 10, /csT/eo 0.02, ^n/eo 20, cJc/eo 100 and n = 0. 



strictly connected. 

A qualitative understanding of this mechanism can 



state with an even number n of electrons. The most 
likely scenario is a spin down electron leaking through 
barrier B into the antidot, since spin up electrons from 
edge A experience a thicker tunneling barrier. The spin 
down electron can now only escape through barrier A, 
and the sequence starts over. The antidot most likely 
lies in the occupied state with 5 = — 1. These processes 
are depicted in Fig. [sja). This mechanism implies a net 
flux of spin down electrons traversing the antidot from 
B to A, leading to a net charge transfer and enabling 
the antidot to operate as a spin filter. Tuning the 
antidot into region //, additional tunneling processes 
are allowed (see Fig. [5|b)). In particular, once a spin 
down electron tunnels into the antidot through 5, a 
s^in up electron can leave the antidot through the same 
barrier. This process, which leads the antidot into a 
state with s = —2, was forbidden in region /, due to 
energy constrains. The antidot is now "trapped" into 
the state with s = —2 until a spin up electron leaks 
through barrier A into the antidot, leading back to the 
state with 5 = — 1, and the sequence starts over. Thus, 
in this regime, a net flow of charge from edge A to edge 
B is created, which leads to a reverse of the current 
with respect to region /. The antidot most likely lies 
in the unoccupied state with 5 = — 2, decreasing the 
mean value of the spin of the antidot. This mechanism 
is repeated when new zero-mode spin excited states 
become allowed. 

In general, the value of Isp depends on asymmetry 
and interactions, as shown in Fig. [Oj where the spin- 
polarized current is plotted for different values of 77 
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and K. In order to explain the detailed form of /gp, 
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Figure 6. (Color online) Plot of hp (units of hVQ/2) as a 
function of V (units of eo/e) for rig = 0.5 with (a-c) K — 1 
and (b-d) K = 0.45. The asymmetry 77 = 1.5 for the upper 
panels (a-b) and 77 = 20 for the lower panels (c-d). In all 
panels /csT/eo = 0.02, En/tQ = 20, cJc/eo = 100 and n = 0. 

plasmon excitations must be taken into account. By 
increasing the bias, more and more collective excitations 
become allowed, eventually making it easier to tunnel 
across barrier A than across even if the former is, 
in principle, thicker. The competition between the 
asymmetry and the presence of collective excitations 
affects the value of the current. From Fig. [6j we note 
that \Isp\ develops a maximum for a certain value of 
bias. If the asymmetry is little, the maximum is at low 
bias; by increasing the asymmetry, the maximum moves 
to higher bias. 

In order to give quantitative analytical results on the 
behavior of /gp^ we have developed a few-states model 
for regions / and //, for the non-interacting case. In this 
regime, the allowed transitions are |n,0) ^ |n + 
(in both regions) and \n + 1, ±1) ^ |n, ±2) (in region // 
only). The tunneling rates (units Tq) are given by 



ri1^,+i(0,0) = ao rfil-i(0,0) = r^ao 

ri1^,+i(0.-2) =^io + «i rL^i^^_i(l,+l) = 7/(ao+xai) 

rL1^^+i(l, -1) = ao + xai rfi^_i(0, +2) = 7?(ao + ai) 
rL1^^+i(l, +1) = xao -1) = TO. (22) 

Here, x = if only states up to \s\ = 1 are involved 
in transport, while x = 1 if states up to \s\ = 2 are. 
From Eq. (14), for K = 1 one has ao = ai = 1. By 



solving the master equation one obtains the stationary 
spin-polarized current and the mean value of the spin of 
the antidot 



-r?+l) 

(772 + 1)2+277(772+^7+1) 



region / 
region //. 



(23) 



-i V _1 

^^+^+l77>l 

2(r7"-l)(r7"+r7+l) , 
(772 + 1)2+277(772+77+1) 



region / 
-2 region //. 



(24) 



As expected, for 77 = 1 no spin-polarized current is 
produced. At 77 7^ 1, spin-polarized current changes 
sign in passing from region / to region //. Thus, the 
system behaves as an electrically controlled source of 
spin currents, whose sign can be reversed by tuning 
the bias. Note that the value of Igp depends on rj: 
for example, for r] < r]* 2.9, |/, 



sp\ 



region / is 

greater than \Isp\ in region //, while the opposite 
happens if 77 > 7^*. This reproduces a maximum of 
\Isp\ moving to higher bias by increasing the asymme- 
try, in agreement with Fig. [6] As already observed, 
the switching of the spin-polarized current is related 
to a decrease of the mean value of the spin of the antidot. 

These results lead to two important conclusions, 
namely that {%) an asymmetric antidot (77 7^ 1) acts as 
a spin filter on the systems, leading to a spin-polarized 
longitudinal current. In addition, (n) the sign of this 
spin-polarized current can be easily switched by tuning 
the bias voltage across the system. Enhancement and/or 
stability in the presence of inter-edge interactions with 
self-consistent screening'^ deserve further investigations. 



V. CONCLUSIONS 

We proposed a device for generating and controlling 
spin-polarized currents in a two-terminal setup, without 
the need of magnetic fields or ferromagnets. The close 
connection between spin and chirality, hallmark of a 
quantum spin-Hall system, is the necessary condition to 
generate this process. 

The spin-filter is realized in a two-terminal antidot 
configuration with asymmetrical tunnel barriers, in 
which helical edge states appear along the edges of the 
bar and around the antidot. By employing a master 
equation approach, appropriate for sequential tunneling 
regime, we have shown how the spin- filtering mechanism 
occurs, and we have related it to a population of the 
high excited antidot spin states. We have discussed the 
role of Coulomb interactions and asymmetry, showing 
how the spin-polarized current can be controlled by 
simple electrical means. 
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